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Abstract 



We reconsider the problem of calculating a general spectral correlation function contain- 
ing an arbitrary number of products and ratios of characteristic polynomials for a, N x N 
random matrix taken from the Gaussian Unitary Ensemble (GUE). Deviating from the stan- 
I dard "supersymmetry" approach, we integrate out Grassmann variables at the early stage 

^ . and circumvent the use of the Hubbard-Stratonovich transformation in the "bosonic" sector. 

1/^ ' The method, suggested recently by one of us is shown to be capable of calculation when 

, reinforced with a generalization of the Itzykson-Zuber integral to a non-compact integration 

■ manifold. We arrive to such a generalisation by discussing the Duistermaat-Heckman localiza- 

tion principle for integrals over non-compact homogeneous Kahler manifolds. In the limit of 
»' , large A'' the asymptotic expression for the correlation function reproduces the result outlined 

' earlier by Andreev and Simons p^ . 



di' 1 Introduction 



Recently there was an outburst of research activity related to investigating the moments and 

correlation functions of characteristic polynomials Zjv(/i) — det ^/ilw — for random N x N 

matrices H of various types. Those studies were motivated by hope to relate statistics of zeroes 
of the Riemann zeta function to that of eigenvalues of large random matrices ||, ||, |3[ |[ |[ ^ , as 
well as by numerous applications of spectral determinants in the theory of quantum chaotic and 
disordered systems §, @, |, |, |l^, |ll|, |l6|, p7|, pl , |l9|, |o), quantum chromodynamics 

I , and relations to interesting combinatorial problems ]23|. 

There are several analytical techniques for dealing with the integer moments (positive or neg- 
ative) of characteristic polynomials. Their applicability varies with the nature of the underlying 
random matrix ensemble. For unitary random matrices ( the so-called "unitary circular ensem- 
ble") one can either relate evaluation of the moments to the Selberg-type integral |jl|, or rely 
upon variants of the character expansion, directly ||2^ or indirectly Another well-studied 

case relates to the ensembles of Hermitian random matrices characterized by unitary invariant 
probability measures. That invariance allows for positive integer moments of the characteristic 
polynomials to be evaluated by methods resorting to orthogonal polynomials The par- 

ticular case of the Gaussian measure can be studied very efficiently by the following procedure. 
First, one represents each of the characteristic polynomials as a Gaussian integral over anticom- 
muting (Grassmann) variables. This allows to average the resulting expressions immediately. At 
the next step one employs the so-called Hubbard-Stratonovich transformation combined with the 
subsequent exploitation of the Itzykson-Zuber -Harich-Chandra integral (24[ p5| : 

f , ^ r / +Nn det |"exp(ix'z/'^)] , , , 

j d^g) exp [t Tr (AgFgt)] = eonst ^^x)\iY) '^^'^'^^^ 

gGU(N)/T 
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Here X, Y are diagonal matrices with eigenvalues x\ correspondingly , where indices k and / take 
values from 1 to A'". A(X), A(F) stand for the Vandermonde determinants and T = U(l) x • • • x U(l) 
is the maximal torus of the group U(N) . A detailed outline of the method can be found, e.g. in 

The evaluation of the negative integer moments of characteristic polynomials turns out to be 
more tricky as care must be taken to avoid divergences and account for presence of poles Jl9| , p2t . 
The standard way goes back to the work by Schafer and Wegner p7[ and relies upon representing 
the (regularized) inverse determinants as the Gaussian integrals over commuting complex variables. 
This method then exploits an extension of the Hubbard- Stratonovich transformation to a non- 
trivial manifold with inherent "hyperbolic" structure. An alternative variant of the method was 
suggested in a recent paper by one of the present authors , referred as [I] henceforth. The latter 
work contains a detailed discussion of the problem as well as many related references. 

A much more general correlation function of characteristic polynomials is one combining pres- 
ence of both positive and negative integer moments. Such a correlation function contains a very 
detailed information about spectra of random matrices and thus it is most important for applica- 
tions in physics. Correlation functions of that type are also interesting for the sake of comparison 
with more refined conjectures on the behavior of the Riemann zeta-function||2^ . 

The standard technique in that case naturally combines Gaussian integrals over commuting 
and anticommuting variables and the subsequent Hubbard-Stratonovich transformation of usual 
and "hyperbolic" nature. The method is known in the literature as the supermatrix (or "super- 
symmetry") approach pioneered by Efetov |p9| in the theory of disordered systems, and taken over 
to random matrices by Verbaarschot and Zirnbauer 

It should work, in principle, for the general case but technically any general calculation be- 
yond the two-point correlation function proved to be extremely difficult. The main problem is 
related to the so-called " anomalous" , or " boundary" terms arising when changing variables in 
the superintegrals. Those anomalies can be traced back to the admixture of nilpotent terms ( 
those containing even number of grassmannian factors) to the usual commuting variables. All such 
anomalous terms were classified by Rothstein[^ in a general form. However, to write down their 
contribution explicitly in a specific parameterization proves to be a very daunting job. The latter 
fact makes the standard supersymmetric calculations usually impractical beyond a few lower-order 
correlation functions. A notable exception is the general many-point correlation function of spec- 
tral densities (see works by Zirnbauer Q and SzaboQ). In the latter case the boundary terms 
do not contribute and the final result is provided solely by the "bulk" integral which is evaluated 
by standard methods. At the same time our main object of interest - the correlator of spectral 
determinants - contains anomalous terms on equal footing with the bulk contribution. 

A modification of the supermatrix method was suggested by Guhr |3^. It relied upon a 
generalization of the Itzykson-Zuber integral to a unitary supergroup and allowed one to go beyond 
the two-point function for the case of Gaussian Unitary Ensemble (GUE). A few years ago Andreev 
and Simons presented an asymptotic formula for a correlation function containing both products 
and ratios of characteristic polynomials for GUE matrices. In their short communicationjl^ they 
claimed that Guhr's method equipped with the further extension of the Itzykson-Zuber integral 
to a psewrfounitary supergroup solved the problem. The authors indicated that the arguments 
behind that generalization were similar in spirit to those by Guhr, but "technically involved". 
They promised to present details of the method "in a longer paper" which, unfortunately, never 
appeared. Let us mention that the role of boundary contributions in the Guhr's method seems to 
be not clearly discussed in the literature. 

In the present paper we show that the method suggested in [I] enables one to calculate the 
general correlation function of integer moments of GUE characteristic polynomials. In fact, we 
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integrate out Grassmann variables at a very early stage thus seriously departing from the general 
spirit of supersymmetry. One of the advantages is that no anomalous terms can ever arise along 
such a route. 

The calculation required, however, the knowledge of an analogue of the Itzykson-Zuber type 
integral over non-compact manifolds related to pseudo-unitary groups. A pseudo-unitary group G 
is defined by the following conditions imposed on the group elements: 

g^Xg = \, V.geG (2) 

Here A is a diagonal matrix with the elements ±1 in a matrix representation. 

In the mathematical literature the integrals of the Itzykson-Zuber type over semi-simple Lie 
groups are interpreted as Fourier transforms over adjoint orbits gXg~^, € G and formulas for such 
integrals are known (see the works of Rossmann |Q, Berline and Vergne [0, Prato and WuQ, 
Paradan . However the derivation of the Itzykson-Zuber type integrals over non-compact 

groups in a manner accessible to physicists was not presented before to the best of our knowledge. 

A particular case of non-compact extension of the Itzykson-Zuber integral was presented re- 
cently in the Appendix C of [I]. The method of calculation, however, relied upon a specific pa- 
rameterization of the integration manifold and seemed hardly applicable to a general non-compact 
case. 

The standard procedure of the derivation of the Itzykson-Zuber type integrals is to employ the 
diffusion equation arguments. This method goes back to the original paper by Itzykson and Zuber 
Q and was also used by Guhr Q, and, apparently by Andreev and Simons Q. We provide such 
derivation for the Itzykson-Zuber type integral over the pseudo-unitary group in our Appendix A. 

In this paper we consider in detail a different approach based on the Duistermaat-Heckman 
localisation principle outlined in the context of random matrices by Zirnbauer |3^] . Indeed, it is 
well known that the Itzykson-Zuber integral is a representative of the family of integrals of the form 
J fi™ exp iH going over a 2m— dimensional phase space. These integrals are semiclassically exact 
provided the Hamiltonian H is "localisable" . In the case of compact phase spaces localisability 
is equivalent to the condition for the phase flow generated by the Hamiltonian H to preserve the 
Riemannian metric of the phase space. In particular, considering the manifold U(N)/T as a phase 
space, the Itzykson-Zuber Hamiltonian 

H{X,Y,g) = Tv{XgYg^) (3) 

appears to be localisable in the sense of Duistermaat and Heckman. 

From the geometrical point of view the underlying phase space U(N)/T of the Itzykson-Zuber 
formula belongs to the family of compact flag manifolds . The flag manifolds are Kahler homoge- 
neous spaces of the type G/H, where G is called a transformation group and H is the centralizer 
of a sub-torus of T in the group G, i.e. H — {g E G\g^^Tog — Tq}, where Tq is a sub-torus of 
T. The compact flag manifolds were considered in details by Picken |Q in the context of the 
Duistermaat-Heckman formula. Picken has obtained an explicit expression for the set of localiz- 
able Hamiltonians on a compact flag manifold. Then the fact that the Itzykson-Zuber Hamiltonian 
is localisable follows as a particular case (G = U(N), H = U(l) x . . . x U(l)). 

The possibility to extend the Duistermaat-Heckman localisation principle to non-compact sym- 
plectic spaces was first discussed by Prato and Wu . The case of the non-compact counterpart 
of the space CP^ was considered by Fujii and Funahashi . However the calculations done in 
[I] suggest to concentrate on more general homogeneous non-compact symplectic manifolds. We 
show that an extension of the Duistermaat-Heckman formula is possible when the homogeneous 
manifold under considerations, G/H is a non-compact Kdhler manifold. The key feature of such 
manifolds enabling one to apply the Duistermaat-Heckman theorem is that their invariant metrics 
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can be chosen to be sign-definite. In turn, the review paper by Bordemann, Forger and Romer 
pO| provides conditions under which a homogeneous manifold turns out to be Kahlerian. The 
conditions can be summarized as follows. Assume that the group G is a connected and semisimple 
and the subgroup H is a compact centralizer of a torus in G. Then 

1. If M = G/H is compact, M is a Kahler manifold. 

2. Let M = G/H be non-compact and L be the maximal compact subgroup of G containing H. 
Then M = G/H is a Kahler manifold if and only if G/L is a Hermitian symmetric space. 

In particular, the manifold U(ni, n2)/T appears to be Kahlerian in this case. Indeed, the maximal 
compact subgroup of U(ni, n2) containing H =TisL = U(ni) x U(n2). Since U(ni, n2)/U(ni) x U(n2) 
is a Hermitian symmetric space of the type AIM (see Helgason p.354) the manifold U(ni, n2)/T 
is a Kahler manifold as follows from the criterions presented above. 

Our general construction is then exploited to derive the Itzykson-Zuber type integral over a 
pseudo-unitary group 

J d^l{g)exp[i^T{XgYg-')]^ 

gGU(ni,n2)/T 

det [exp(zx'y'=)] |i<,_fc<„^ det [expiix^yf")] |„^+i<,,fc<„^+„^ 

" 'A(x)A(y) " 

It then can be seen that the remarkable formula of Harish-Chandra remains valid also for the 
non-compact group U(ni, 02). Indeed, the righthand part of Eq.(^) can be rewritten as 



a>0 a>0 

where W is the Weyl group corresponding to U(ni,n2), i.e. W = Sp^ x Spj. The a{X) is a root 
corresponding to a Cartan subalgebra element X, w{Y) — wYw^^ and \w\ denotes the parity of 
w. 

The structure of the paper is as follows. In sections 2 and 3 we provide a necessary back- 
ground information on Kahler geometry and discuss the Duistermaat-Heckman localisation on 
non-compact Kahler manifolds. Then in section 4 we give an account of (a refined version of) the 
method suggested in [I] for calculating the general correlation function of characteristic polyno- 
mials for the GUE matrices. In section 5 we analyse the derived matrix integral representation. 
For this purpose we use the formulae Eqs.(^ and (^) and evaluate the remaining integrals by the 
saddle-point method in the limit N 00. The open questions are summarized in the Conclusions. 
Technical details are presented in the appendices. 



2 Basic properties of Kahler manifolds 

In this section we show how the Duistermaat-Heckman localisation principle can be refor- 
mulated for the Kahlerian dynamical systems. These are dynamical systems whose phase spaces 
are (simply connected) homogeneous Kahler manifolds ( for definitions and basic properties of 
homogeneous Kahler manifolds see, for example, Kobayashi and Nomizu ||49[|). 

Throughout the paper we use a complex parameterization on flag manifolds which is introduced 
following Borel's method Q . A detailed exposition of the method can be found in the papers by 
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Bar-Moshe and Marinov 1 47 , |Q . Below we provide a concise description of the main features of 
the underlying structures. 

Given a semi-simple Lie group G we introduce the canonical Cartan-Weil basis for the corre- 



sponding complex semi-simple Lie algebra g: {to} = {hj,e±q}. Here a = 1,2, 



dim g; 



j = 1,2, . . . ,r — rank g, and {q} G Ag are the positive roots of the Lie algebra g. The complex 
parameters which are introduced in the flag manifold G/H correspond to the positive roots of the 
Lie algebra g. When the subgroup H of the group G is a maximal torus, H = T, any element g{z, z) 
of the coset space G/T has the following decomposition: 



u(z) — exp 



j{z,z) = u{z)p{z,z) 

( 



, p[z, z) = exp 



E exp ^W{z,z)h^ \ (5) 



where we used the bar to denote the complex conjugation. 

To find an explicit expression for y'^ and W as functions of the complex coordinates z'^ ,z'^ = z'i , 
one has to exploit the condition g^iz, z) — g~^{z, z) when the group G is unitary, or the constraint 

for the case of a pseudo-unitary group. The element g{z, z) = u{z)p(z, z) constructed in this 
way represents a point with the coordinates (z"^, z^) on the flag manifold G/H. To provide a reader 
with a simple but informative example we show in the Appendix B an application of the general 
principles outlined in this section for the case of the compact manifold U(2)/U(l) x U(l) and its 
non-compact counterpart U(l,l)/U(l) x U(l). 

Let us note that a complex parameterization introduced above can be looked at as a convenient 
way of describing the action of the transformation group G on its flag manifold G/H. Indeed, for 
any g g G the (unique) decomposition gu{z) = u{gz)p{z, g) allows one to flnd gz and thus to 
determine the (holoniorphic) action of the element g of the transformation group G on a point of 
the flag manifold G/H. 

The homogeneous Kahler manifold M comes with the Kahler potential K{z,z), which is a 
scalar function deflned on any open neighborhood of M with local complex coordinates z", {a = 
1,2, . . . ,m ^ dimcM). When a group G acts holomorphically on M, 



gz,yg e G 



the Kahler potential K{z, z) is transformed as 



K{z, z) -> K{gz,gz) = K{z, z) + -f 5 G G 



(6) 



(7) 



where '^giz) is a holomorphic function of z. Once the Kahler potential is provided, the exact (1,1) 
differential form on the Kahler manifold is introduced as follows: |^ 



cu^a {z, z) dz°' A dz'^ , uj^p {z, z) = ——— dadaK{z, z 



2m 



(8) 



where the factor — l/(27ri) is chosen for convenience. Equations Eqs.(^), (|^) show that the (1,1) 
form D. is invariant under the holomorphic action of the group G on the homogeneous Kahler 
manifold M. When the phase space (Kahler manifold M) and the (1,1) form are specified, the 
classical mechanics is defined by the Poisson brackets for any two smooth functions Fi{z, z) and 
F2{z,z): 



{F^{z, z),F2{z, z)}p5. = w'^'^iz, z) {d^Fiiz, z)dpF2{z, z) - d^F^iz, z)dpF^{z, z)) 



(9) 



^Throughout the paper we assume tensor notations, i.e. summation over repeating indices 
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where the antisymmetric field w"^{z,z) is inverse to w^p{z,z), i.e. 

w°''^{z,z)w^,p{z,z) ^ 6"^, w''^{z,z)w^a{z,z) ^ (10) 

Provided that a Hamiltonian H{z, z) is given , the equations of motion on the Kahler homogeneous 
manifold M can be written in terms of the Poisson brackets, 

dF{z,z)/dt^ {F,H}pg (11) 

From Eqs.(0)-(|ll|) we then obtain a time evolution of the complex coordinates: 

z'' = w''^{z,z)d0H, z^ = -w''^{z,z)da,H (12) 

In the formulation of the Duistermaat-Heckman localisation principle one uses an important 
notion of the Hamiltonian phase flow defined on the homogeneous Kahler manifold M as a one- 
parameter group g{t) of diffeomorphisms M ^ M (see, for example, Arnold [^): 

d/dt\t=o (git) ■ 2") = w°''^{z, z)d0H, d/dt\t=o [g{t) ■ z^) - -w^^iz, z)daH (13) 

From the above definition and Eqs.(|7|), (^ it follows that the Kahler metric remains invariant 
under the holomorphic phase flow. 

Let us now represent an arbitrary element of the Lie group G acting holomorphically on the 
Kahler manifold M by the set of cartesian coordinates: 

g(0 = exp (^VJ , a = l,2,...,n = dimg (14) 

where are basis elements of the Lie algebra g of the group G satisfying the commutation relations 

[Ta,n]= fabTc, (15) 

with /^^ being the structure constants of the Lie algebra g. 

The Lie group G is itself a homogeneous space so the left action of the group on itself induces 
vector fields: 

Ta - D,(e) - LliOdi (16) 

where denotes a partial derivative with respect to the parameter 

The holomorphic action of the Lie group G on the Kahler manifold M z G M ^ gz G M, G G 
induces the vector fields of the form 

Va{z,z)^K2{z)d^ + K^{z)da, (17) 

where the fields k"(z) are expressed by the induced vector fields on the group: 

<{z) = Da{0{9m"k=o (18) 

The conjugate fields k"{z) are defined similarly. Once the Va(z,z) are induced vector fields they 
also satisfy the commutation relations (|l^) of the Lie algebra g. 

Introduce now the linear operators ia acting on the differential form Q as: 

ia^ = w^p{z, z)K"{z)dz'^ - w^^{z, z)n^{z)dz" (19) 
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The linear operators ia are related to the Lie derivatives Lv„ corresponding to the vector fields 

Lv„ = dia + iad (20) 

where d stands for the usual external derivative operator acting on forms. The invariance of the 
(1,1) Kahler form fl, Eq. (||) under the holomorphic action of the Lie group G can be represented 
by the condition 

Lv„rj = 0, Va = 1,2, . . . ,71 = dim g (21) 
Since the Kahler form is closed: = the condition Eq.(^l|) can be rewritten as 

dia^ = (22) 

As the Kahler manifold M is assumed to be simply connected, and in simply connected spaces all 
closed forms are exact: dfl = ^ D, ~ diu, the equation above implies the existence of n = dim g 
functions Ta{z,z) (unique up to an additive constant) such that 

iM = {2TTi)-^dTa (23) 



In terms of the local complex coordinates on the Kahler manifold M differential equations ( p3D 
acquire the form: 

daTa{z, z) = -27riKf (z)w„^(z, z), dfjTa{z, z) = 2TriK"{z)w^0{z, z) (24) 
An immediate consequence of these equations is 

w"'^(z, z)dMz, z) = -2m4{z). uj^f^iz, z)d-pTa{z, z) = 2^i<(z) (25) 

which, in turn, implies the relations: 

dp {uj^^'iz, z)d,Ta{z, z)) = 0, a„ (w''^(z, z)d^Ta{z, z)) = (26) 

The functions ^^(z, z) are called equivariant momentum maps, and we will see that they play 
an important role in providing the relation between the Kahler geometry and the Duistermaat- 
Heckman localisation principle. The explicit construction of the equivariant momentum maps 
Ta(z, z) for the compact flag manifolds was performed by Bar-Moshe and Marinovp?] and will 
be discussed later on in the present paper. 

2.1 Duistermaat-Heckman localisation and equivariant momentum maps 

Now we will obtain the conditions under which a Hamiltonian i?(z, z) on the Kahler manifold 
can be localisable in the sense of the Duistermaat-Heckman theorem. In this section we adopt 
the method developed previously by Bismut ^\ , Witten|52j, Zirnbauerpof and apply it for the 
particular case of the Kahler manifolds. 
We consider the integral 

''n'"exp(iF(z,z)) (27) 



where M is a Kahler manifold on which the complex parameterization is introduced. When M is a 
flag manifold it has a complex parameterization as discussed above, and the (1,1) form O is given 
by Eq.(|^), with m being the complex dimension of the manifold M, dim^M = m. Let us introduce 



7 



2m anticommuting variables that are counterparts to commuting complex coordinates 

(z", z") of the manifold M. Integral ( p7| ) can be rewritten as that with a flat integration measure: 

/m m 
Jl dz"dz" Yl dfdfexpS{z,z,^,^) (28) 

a,a=l /3j=l 

where 

^(z, z, = z) - w^p{z, z)Ce (29) 

We will refer to the expression in the exponent as to the "action" depending on the commuting 
complex coordinates (z", z") of the manifold M and on the anticommuting variables (^"I'C")- 
Let us introduce a first-order differential operator D defined by the formula 

D = Cda + Cdo. - tw'^^iz, z) (d^H{z, z)dj - d^H{z, z)9«) (30) 

Using explicit expressions given above one can verify that such a differential operator annihilates 
the action S, i.e. DS = 0. Next step is to construct a function A(z,z, ^,^) on the extended space 
which is annihilated by a repeated action of the operator D, i.e. D^X = 0. When such a function 
exists and integral ( |2^ ) converges it is possible to deform integral (|2^) as follows, 

/rn m 
Yl dz^'dz" Y[ d£,'^d£,^exp{S{z,z,^,^)+WX{z,z,£_,^)) (31) 

a,a=l /3,/3=l 

where t is an arbitrary parameter. Indeed, one can expand the integrand in a series with respect 
to the parameter t and use the integration by parts together with the properties of the differential 
operator D {DS = 0, D^A = 0) to verify that the integral It does not depend on the parameter 
t, i.e It = I- Following the general procedure described by Zirnbauer[^ we make the following 
choice for the function A(z, z, ^, ^): 

A(z, z, ^, - t {d„H{z, z)r - d^H{z, -z)e) (32) 

The action of the first-order differential operator D on this function gives 

D\{z, z, = -2 (w^\z, -z)d^H{z, z)dpH{z, z) + id^pH{z, z)^^'^) (33) 

A repeated action of the operator D on the function A gives the following expression: 

i?2A(z,z,C,0 = -2r5a {w^''{z,z)d^H{z,z))+2edc. {w''''{z,z)d,H{z,z)) (34) 

It immediately follows that the function A(z, z, ^, ^) satisfies D^X = if and only if the Hamiltonian 
i?(z, z) satisfies the following conditions: 

da (w'^'^iz, z)d^,H{z, z)) = 0, dp {w''''{z, z)d^H{z, z)) = (35) 

The property of the Kahler cosets that ensures localisation is that one can always choose on them 
a sign-definite Riemann metric (see, for example, Kobayashi and Nomizu p9| and a review article 
by Bordemann, Forger and Romer |40|). When the metric is positive definite the numerical part 
of the expression for DX proportional to 

{dH)^ = u;"^(z, z)daH{z, z)di^H{z, z) (36) 
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is positive definite as well. It then follows that the limit t ^ oo localises the integral It in Eq. (^Tj) 
on the critical set dH = 0. In turn, it implies that the original integral / is localized on the critical 
set of the Hamiltonian H as well. For the negative definite case one can just set t —oo with the 
same result. 

As is seen from the definition of the Hamiltonian phase flow on a Kahler manifold, Eq.(|l3|), con- 
ditions Eq.(|35|) mean that the phase flow generated by the localizable Hamiltonian is holomorphic. 
Such phase flow preserves the Kahler metric, as follows from definition (|l3|) . 

In particular, consider a Hamiltonian on the Kahler manifold M which can be represented as a 
linear combination of the momentum maps Ta{z, z) satisfying the Eqs.(|2^), 

n 

= ^c,T„(z,z) (37) 



Then the relation Eq.(|25|) ensures that such Hamiltonians conform to the conditions (|35[ ). Hence 
they are localisable in the sense of the Duistermaat-Heckman principle provided the integral ( p7| ) 
converges and {dHY defined by equation is sign-definite. 



2.2 Localisable Hamiltonians on flag manifolds with unitary transfor- 
mation group 

Let us consider first the case of the flag manifold G/H, with G being a unitary transformation group, 
Gt — G^-"-. In any matrix representation there exist projection matrices r^j which correspond to the 
elements hj of the Cartan subalgebra of the Lie algebra g. The projection matrices are defined by 
the following set of equations 

V] = = Vj Vjhk = hkVj Vj, /c = 1, . . . , r 

Vje^qVj = e^qVj VjeqVj = Vj^q (38) 

Here the hat stands for the matrix representation. When the Cartan subalgebra elements are 
represented by iV x diagonal matrices, the projection matrices are also diagonal, and any N x N 
diagonal matrix is a linear combination of the projection matrices and the unit N x N matrix, Ijy- 

The projection matrices were introduced by Bando, Kuratomo, Maskawa and Uehara 
and Itoh, Kugo and Kunitomo |Q to construct explicit formulae for the Kahler potentials. In 
particular it was found that in the case of a flag manifold with a unitary transformation group 
G the most general Kahler potential is a linear combination of scalar functions (the fundamental 
Kahler potentials). Each fundamental Kahler potential corresponds to a basis element of the 
Cartan subalgebra of the Lie algebra of the group G. The explicit expression for the fundamental 
Kahler potential corresponding to the basis element hi of the Cartan subalgebra is given by 

K,{z,z) = In det {riiU^z)u{z)rit + I - Tj,) , m(z) = exp(z'^eg) (39) 

This formula was used by Bar-Moshe and Marinov ^ to derive the equivariant momentum 
maps in terms of the local complex coordinates on a Kahler manifold. The expression obtained by 
Bar-Moshe and Marinov is 

r 

Ta{z,z) ^ -TT{p{z,z)fa) , p(z,z) = ^liPi{z,z) (40) 

1=1 

where li are arbitrary constant coefficients and the matrices pi{z,z) are given by the following 
equation: 

Pi{z,z) ^ u{z)rji{rjiU^ {z)u{z)r]i + I - rji) ^ rjiU^z), u{z) = exp {z'' iq) (41) 
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The Hamiltonians on Kahler manifolds which can be represented as a hnear combination of the 
momentum maps are locahsable. For the corresponding integrals the Duistermaat-Heckman for- 
mula is applicable as we have seen in the previous section. Using the expressions Eq.(^ for the 
momentum maps we thus obtain the following formula for the localisable Hamiltonians: 

n r 

H{z, z)=Y,J2 "«Tr {p^{z, z)Ta) , (42) 

a— 1 i—1 

where Vai are arbitrary constant coefhcients. The matrices pi{z,z) are transformed under the 
action of the group G on the manifold M as 

Pi{z,z) ^ p,{gz,gz) = gp^(z,z)g\ = g^^ (43) 

The above transformation law for the matrices pi{z,z) can be verified using the decomposition 
gu{z) = u{gz)p{z, g) which determines the action of the transformation group G on its flag manifold 
(see Appendix C for explicit calculation). Note that pi(0, 0) = rji and the projection matrix pi{z, z) 
at the point with the complex coordinates (z", z") can be written as 

p^{z,z) = g{z,z)i},g-'^{z,z). (44) 

Here the group element g{z,z) e G represents a point of the coset space G/H parameterized as 
described earlier in this section. 

A short inspection of expressions Eq.(^2[) and Eq.(|4^) for the particular case G — li{N) and 
H = T makes it clear that the localizable Hamiltonian H{z, z) is the same as that entering the 
Itzykson-Zuber integral, see Eq. (||). 

The construction described above can be taken over to the case of pseudo-unitary transforma- 
tion group without much modification. 



3 Duistermaat-Heckman localisation principle for manifolds 
with pseudo-unitary transformation group 

We consider the flag manifolds G /T where the transformation group G is a pseudo- unitary Lie group 
with its elements satisfying the condition (||). Localisable Hamiltonians will be expressed as linear 
combinations of momentum maps similar to the case of the unitary transformation groups. The 
only new element is the presence of the matrix A. We begin with the formula for the fundamental 
Kahler potentials (compare with Eq. ([39|)): 

Ki{z, z) = In det [riiu\z)Xu{z)rii + 1 - rji) , u(z) = exp(z^eg) (45) 

Once the fundamental Kahler potentials are known it is possible to find the equivariant momentum 
maps in terms of the local complex coordinates on the given flag manifold. We have found that 
the equivariant maps have the same form as in the unitary case (see Eq.(pO[)), but the projection 
matrices Pi{z, z) turned out to be slightly different and are given by 

Pi{z, z) = u{z)rii {j-jiU^ {z)\u{z)r]i -fl - rji) ^ rj^u^ {z)\ (46) 

We note that in the pseudounitary case the element u^{z) in the above expression comes together 
with the matrix A in the same way as it does in the formula Eq.(^5|) for the fundamental Kahler 
potentials. Respectively, the Hamiltonians given by the formula Eq.(^2[) with pi{z,z) specified by 
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Eq.([45|) are localisable provided that the integral Eq.(|2j) converges. The transformation law for 
the matrices Pi{z, z) turns out to be of the same form as one for the unitary case: 

p^{z,z) ^ p,{gz,-gz) = gpi{z,z)g^^, g^ Xg = X (47) 

Then the same argumentation makes it evident that the natural counterpart of the Itzykson-Zuber 
Hamiltonian is 

Hig,X,Y)=Tr{g{z,z)Xg~\z,z)Y), .g(z, z) G U(ni, 7^2)/T (48) 

This Hamiltonian is localisable provided that the diagonal matrices X, Y are chosen properly as 
to ensure that the integral Eq.(^) converges. 

In what follows we apply the Duistermaat-Heckman theorem to the integral Eq.(|^) with the 
Hamiltonian given by the expression Eq.(p8[) and evaluate it by using the method of stationary 
phase. We conventionally refer to this procedure as to the " semiclassical approximation" , but the 
Duistermaat-Heckman localisation principle ensures that such an approximation yields the exact 
result. 

The elements g of the coset space U(ni, n2)/T satisfy the condition Eq.(||) where G — U(ni, n2), 
and the matrix A is given by 

'-{'o -I) <«> 

For our calculation we make use of the same complex parameterization for the manifold G — 
U(ni, n2)/T as that described in section 2. In the integral Eq.(^) the matrices X,Y are diagonal, 
with the complex variables Xi, . . . , Xni+n2 ^-nd j/i, . . . , yni+n2 ■ 

Let us find now the set of solutions for the equation dH ~ (saddle points). We get: 

dH = Tr {Sg [Y,g-'Xg]) =0, Sg = g-'dg (50) 

Note that the expression Tr(x, y), with x e g, y € g is a quadratic form on the Lie algebra g. 
Then from the condition Tr(x, y) = follows that x = or y = 0. We see that the equation ( |50| ) is 
equivalent to the condition of the Lie algebra elements Y and g^^Xg commuting with each other, 
i.e. 

[Y,g-'Xg]^0 (51) 

Once the matrices X, Y are diagonal all the saddle points on the homogeneous manifold 
U(ni, n2)/T belong to the permutation group Snj+n2 (which is the permutation group of the diago- 
nal entries of X and Y). Denote these saddle points by P. The saddle points P should belong to 
the coset space U(ni, n2)/H so they are elements of the pseudo-unitary group U(ni, 02) and satisfy 
the following constraint: 

P^XP = X, PeSni+n, (52) 

It is clear that only permutation matrices that do not mix the elements -1-1 and —1 of the matrix 
A could satisfy the condition Eq.(^2|). Thus the relevant saddle points are all elements of the 
permutation group Spj x Spj, where Sm stands for the permutation group of the top rii diagonal 
elements of the matrix A, and Snj is the permutation group of the rest n2 diagonal entries of A. 
Let us recall that in the case of unitary coset space the relevant saddle points were all possible 
elements of Sm+nj. We conclude that the condition ensuring the saddle points to be elements of 
the pseudo-unitary group reduces the initial symmetry group Sm+nj of saddle points down to the 
symmetry group of the lowest order, x ■ 
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In order to evaluate the integral Eq. (Q) by the stationary phase method we rewrite it as a sum 
over contributions from neighborhoods of the saddle points 

f d/i(g)exp(zi/(5,X,y))= ^ ( dti[gj,)e^^[iH{gp,X,Y)) (53) 



geU(ni,n2)/T 



PeSnj xs„. 



where gp = g ■ P. Each integral in the sum above corresponds to one saddle point P and should 
be taken over a neighborhood of the point P on the manifold U(ni,n2)/T. Consider the shift 
gp — g - P as a change of coordinates. Then it suffices to take the shift element g to be close to the 
unit element of the group. This allows one to explore the neighborhood of the chosen saddle point 
solution in line with the spirit of the Duistermaat-Heckman principle. Introducing the complex 
parameterization g = g{z, z) we therefore take the complex parameters of the shift element g to 
be close to zero. We note that the coset measure is invariant under the group shifts, 

dfj.{gp) = dfi{g ■ P) = dfi{g) (54) 

enabling us to use the same flat measure at each integral in the sum in the expression Eq.(^3|). 
Correspondingly, the Hamiltonian H{gp, X, Y) is expanded in the vicinity of the relevant saddle 
point P separately at each of the integrals in the right hand side of Eq.(p3|). 

For each g G U(ni,n2)/T represented by a point from the coset space we introduce the local 
complex parameterization in the neighborhood of unity, g = /. We have: 

5(0,0) = /, giz,z) = l + g^'\z,z)+g^^\z,z) + 0{\z\^) (55) 

(Here (jr(-^^(z,z) = 0(|2:|) and g^^'>{z,z) = 0(|zp)). The condition of pseudounitarity of the trans- 
formation group underlying the given coset space yields: 

g-\z, -z)=I + \ (g^^\z, -z)))\ + A.g(2) -,)x + O {\z\^) (56) 

Then it is straightforward to find the Lie algebra element g(^'(z,z) up to the lowest order in 
z, z using the Cartan-Weyl basis of the Lie algebra g: 

5(i)(z,z) = Z + iV + Q 
Z = z«eq, N^m{z,z)e-q, Q=Q^{z,z)h, (57) 

The Lie algebra elements TV and Q turn out to be of the first order in \z\. Insert Eqs.(|55|), ( ^6| ) 
with g'''^\z,z) given by Eq.(|57|) to the pseudo-unitary condition Eq.(|^). We neglect the terms of 
the second order in |z| and find the Lie algebra elements N and Q: 

N = -AZU, Q = (58) 

We see that the first-order Lie algebra elements in terms of the local complex coordinates are given 
by: 

c,(i)(z,z) = AZU = Z- Zj; (59) 

2 

Z = ^ z'^eq = ^ %e(y) (60) 

g=l l<2<j <ni +Ti2 
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where the standard basis elements e(y) of the Lie algebra are defined in the matrix representation 
by the formula: 

{^i^3)}ki^^^kSji (61) 

Now we should expand the Hamiltonian 

H [gpiz, z)) = tr {Xgp{z, z)Ygp\z, z)) (62) 

in the vicinity of the saddle point P. For this purpose we insert Eqs.(p5|), ( ^ ) into the Hamiltonian 
and obtain 

H (5p(z, z)) = + (z, z) + (z, z) + O (|z|3) 
= tr {XYp) 

(1) (z, z) = tr (^XYp (.gW (z, z)) V 5^'^ 5)' 

(2) (z, z) = tr (^Xg^i) (z, z)rp (g^ (z, z)) + tr (^XYp (gf ^ (z, z)) V y^^) g)^ ^ (63) 

where we have used the following notations: 

(g«)t^A(gW)tA, (^f )t ^ A(g(2))tA, Fp ^ prp^i (64) 

The explicit forms for H'^^^ and iJ*^^^ in terms of the local complex coordinates can be found from 
the pseudo-unitary condition (||). We find that ij'^-* = and iJ*^^-* can be written as 

if (2) (z, z) - tr (^Xgd) (z, z)yp (g « (z, z))^^ - tr (^XFp (g « (z, z)) \9(i) (z, z)^ ^ (65) 

We further observe that taking into account terms of first and second order with respect to |z| re- 
sults in an expression for the Hamiltonian H{gp(z, z)) which depends only on the element g^^-'(z, z). 
Once the expression for g^^\z,z) is found (Eq.(|59|)) we obtain from Eqs.(|63|) and (|65|): 

H{gp{z, z)) = tr (XYp) + tr (zj^X [ZYp] + ZX zIYp ) + O {\z\^) (66) 

ni+n2 ni+n2 

Inserting X = ^ a;*e(ii), Yp = ^ yp^iu) Si,iid Z — ^ Zije(^ij-j to the above equation 

2—1 i— 1 l<i<j<ni+n2 

and calculating the traces yields the Hamiltonian H{gp{z, z)) up to the second order terms in |z|: 

H{gp{z, z)) = x'y'p - ^ (a;' - x^){y'p - y^p)z,jZ,j 

+ J2 (x' -x^Wp-y'p)z,jZ,j- J2 (a;'-a^-'')(yp-J/p)%%+0(|z|3) (67) 

l<z<ni <j <ni +712 ni + l<i<j<ni+n2 

We use the above expansion of the Hamiltonian H(gp{z, z)) in the vicinity of the saddle point P 
to obtain the semiclassical approximation for the integral Eq.(^): 

/= ^ expli^x^y^p] I Y\ ^f!^!lexp[-i ^ {x' - x^){y'p - y^p)z,jZi 

P£S„^xS^2 \ j = l J l<i<j<ni ^ ' \ l<i<j<ni 
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ni + l<'i<j<ni +n2 \ ni + l<i<j<ni+n2 



l<i<ni <j'<ni+n2 \ l<'i<ni <j<ni +n2 

Calculating the Gaussian integrals and noting that for the Vandermonde determinantal factors 
A(F) — ni<j(!/i ~ Vj) O'^^ A (Yp) — {—Y A (y), with e being equal to unity (zero) for even 
(odd) permutations, respectively, after straightforward manipulations we obtain the final formula: 

7- _ (- 1 (exp(^a;^7/')) |i<fc,;<ni det (exp(za;^y')) |ni+i<fc,i<«i+«2 

^ ^ ^ ' A (X) A (r) ^ ^ 

4 Correlation function of characteristic polynomials: gen- 
eral formalism 

In the present section we derive an integral representation for the correlation function of character- 
istic polynomials of the GUE matrices suitable for further investigation in the limit of large matrix 
dimensions. Our method is a refined version of that introduced in [I]. 

Let H he N X N random Hermitian matrix H = W which is characterized by the standard 
(GUE) joint probability density: 

V{H) = Cwexp-yTri?^, Cn = (27r) (70) 

with respect to the measure dH — J^^j^ dHu Ili<j dHijdHij , where as usual dzdz = 2dRezdlm.z . 

Regularizing the characteristic polynomial Z]si{y) = det {^^\n — by considering the spectral 
parameter such that Im/x ^ we are interested in calculating the following generating function: 

\ ri/^l ZN{f^iB)ZN{^^2B) I GUE 

where we denote by (...) the expectation value with respect to the distribution Eq.([70|) and 

jiB = diag(^^^j, ^) , jiF = diag(^^^j, ...,^^'^''^) 

and Im(/z[']j,-/x^']j) > 0. 

The generating function is obviously an analytic one with respect to the complex variables 

{n^ip , IJ'^2f) ■ turns out to be technically convenient to change: /i^^ —il^i^p ' t^2F ~^ ~'^^^2F 
when performing the ensemble averaging, and restore the original generating function by a simple 
analytical continuation. 

To calculate the average we first use the standard " supersymmetrisation" procedure and rep- 
resent each of the characteristic polynomials in the denominantor as the Gaussian integrals: 

[Z^(m)]-i = ^^Id^Sexpjill^^StS-StiJs)! (72) 
Yrf^Sexpji^ (^iS^S-Tr ® St ) | 



(47ri) 



N 
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where we introduced a complex iV— dimensional vector S = (si, sjv)"^ (here T stands for the 
vector transposition) so that (PS — HiLi dsidsi. The sign ± in the exponential is coordinated 
with the sign of Im/i as to ensure convergency of the integral. 

For the characteristic polynomials in the numerator we use the representation in terms of the 
Gaussian integrals over anticommuting (Grassmannian) A^— component vectors x, X^i see e.g. [^ . 
Taking the product of all the integrals, the generating function can be written down in the following 
form: 

/np TLF f 1 / 



d S/,1 / d S;,2 exp ■ 



1=1 



X ( exp 



. 1=1 



(73) 



GUE 



The ensemble average is then easy to perform via the identity: 



-iTr[Hi] 



-5^Tr[A^] 



GUE 



and after straightforward manipulations one arrives at: 



ICN{fib,P-f) oc 



/ Yldxk.idxli / YldxkadxU^^P \ , (^i^^^^Li^'^a + Ai2Fxi^2Xfe,2) + ^Tr [Q 

k k V k 

{, riF 2 / Kb \ 

-i]^ E E Ap E ® ^1 - ® St 2 xk.p 
fc=ip=i \;=i / 



(74) 



X exp 

I. k=ip=i \i=i 

where we introduced the following 2nF x 2nF Hermitian matrices Qp and 2nB x 2nB Hermitian 
matrices Qb- 



with entries 



~ ^li,pi^l2,P2 



^(9l«2) 



and used the notation L — diag(l„^, — 1„^) 

Now we employ the Hubbard-Stratonovich identity: 



fclfc2 



— Xki,qiXk2,q2 



J dQ^e.^ [-f TrQl^ - ^ixli^xUM ( ) 
where the integration goes over the manifold oi2np'x2np Hermitian matrices Qf = 

\<1f If 

with the symmetry structure inherited from that of Qf and wc introduced the shorthand notation: 



(75) 

.(11) .(12) N 

qp rp 

421) .(22) 

9f If 
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Exploiting the above identity allows one to perform the gaussian Grassmannian integral explic- 
itly and to bring the expression to the form: 



1=1 I 1=1 



a (0 c 



X e 



where 



V «l,2nFliV a2,2nFljV 

,(fe) 



Ol,2nFlAr \ 



(76) 



(77) 



We have used notations afei,fe2 = f-^p '^kiM ^ i<lF)kik2 j for k — 1,2, np and introduced the N xN 

matrix B = E?=i [S/,i ^ S/^^ - S/,2 <E> S,^ 

To bring the determinant of the matrix in the above expession to the form suitable for 
further manipulations we consider the case N >2nB and further introduce the N x N matrix 



N 



= {m. 



61,62, . .•,eN-2nB 



where the N-component orthonormal vectors 61,82, ■••,eN-2nB chosen to form a basis of the 
orthogonal complement to the hnear span of vectors Si_i, S„^,i, 81^2, S„^,2 (without restrict- 
ing generality we can consider the latter vectors to be linear independent). Correspondingly, the 
matrix M is chosen to have 2nB columns which arc just the vectors S;^i and Si^2- Simple calculation 
shows that 

M'^M = Qb and M^BM = ^QsLQb (78) 



Then we can write: 



det 



= det 



/ Mj, 
Mi 



N 



^fei ,k2 



Qi 



( 



det 



N-2nE 



ai,i 

02,1 



Ml / 



LjV-2nB 

ai,2 
a2,2 



X det(AB) X det 

QbLQb 



2N 

( Mn 
Mn 



\ 



(79) 



Mn / 



1 

2iV 

ai,2nF \ 
«l,2raF 



'-N-2nB 



1 A 



det ( aki,k2QB - Ski,k2-^QBLQB 



\ 0'l,2nF 0'2,2nF ■■■ 0,2nF,2nF / 

Now it is easy to see that the determinant factor detAg is given by 



det 



Af - Qf 



where we introduced the notation qj^p for (real) eigenvalues of the (Hermitian) matrix Q'"p' = 
fj-F - Qf- 



N-2r. 



2nF 



JJdet 



fc=i 



1 ^ 



(m) _ 
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Next step is to deal with the integrals over S/^i,S;^2- For this we observe that the integrand 
depends on those variables only via the matrix Qb and employ the following 
Theorem I 

Consider a function F(Si, ...,Sm) of A^-component complex vectors S/ \ <l <m such that 



/ d^Si... / d2s„|F(Si,...,S„,)| <oo 



(80) 



Suppose further that the function F depends only on w? scalar products Sj^S/^ 1 < ^1,^2 < so that 
it can be rewritten as a function T{Qm) of m x m Hermitian matrix Qm'- 



I SjSi S{S2 
S^Si 

V st„Si ... 



Then for > to 
where 



rf^S,„F(Si, Sm) — Cat ,,! 

(27r)^"~ 



828™ 
dQm (detQ 



N-r 



HQ^n) (81) 



c 



and the integration in the right-hand side of Eq.(81) goes over the manifold of Hermitian positive 
definite to x to matrices Qm- 

In fact, the formula Eq.(^) was already implicitly used in [I]. In that paper it was justified by 
heuristic arguments employing the Fourier transform of the function J- and subsequent exploitation 
of a matrix integral close to that considered by Ingham and Siegel[|2|. A proof of the theorem is 
given in the Appendix D of the present papeij^. 

In our particular case to = 2nB and the role of Qm is played by Qb- The convergency 
condition Eq.(|80|) in Eq.([76|) is ensured by imaginary parts of the spectral parameters fj.. Applying 
the theorem, we get: 



N-2n 



K,N{fiB, ftp) J dQp (det Q^^^ 
/ dQBdctQ^-2"-e-5^Tr(Q«L)V|Tr[/i«QBL] -Qdet 



(82) 



Let us now replace the Hermitian matrix (3_f by Q^jf' as the integration manifold and further 



change: Qb 2NQb- Omitting both tilde and (fi) symbols henceforth we arrive at: 



/Cjv(Ab, Af) oc / dQp {detQpf'^"^ e-*Tr[A. 



dQs idetQef''"'^ Tr(Q«L)V.ivTr[/,,Q,L] 



Qb>0 



Qi 



(83) 



^Afternomplction of our work we learned that an equivalent formula was used earlier by David, Duplantier and 
Guitter feal in quite a different context 
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The integral representation Eq.(g3|) is our main result for the present section. It is valid for 
any parameters N,nB,F,i^B,F provided N > 2nB- The form of the integrals is clearly suggestive 
of treating them by the saddle-point method in the limit of large N . The details of the procedure 
are described in the following section. 



5 Large behavior of the correlation function 

To perform the saddle-point evaluation of the integrals we first have to expose those degrees of free- 
dom which are amenable to such a treatment. It is immediately evident that for the matrix Qp the 
relevant variables are real eigenvalues —oo<qk<oo^l<k<2np. Accordingly, we write Qf = 
UqpU^ where U € U{2nF) is 2nF x 2nF unitary matrix, and qp = diag(gi, 92nir)- The integra- 
tion measure in those variables is known to be written by dQp oc A^{(ji?}(i/i(L/)(i(7F, with d^iilJ) 
being the corresponding Haar's measure on the group U(2nF) and A{qF} — nfci<fc2 ('^'=1 ^ '^'^'2) 
standing for the Vandermonde determinant (see e.g. [is)). 

The only term in the integrand of Eq. ( ^3|) which depends on the unitary matrix U is obviously 
the exponential expiVTr (/tFC^^FC^^)- We immediately see that the corresponding integral over 
the unitary group is exactly that by Itzyson-Zuber-Harish-Chandra, Eq.(|l|). This yields: 

ICN{liB,M « TT^e"*^"""'' / dqpAiq} {dctqpf-'''" det Um^^«'=2 1'"" (84) 
A{^f| J L Jfci,/C2=i 



A{Af} 

X e-f Tr[5H^ / dQs (detQ^)^"'"" e-^Tr(Q.L)V. A^Tr[A.Q«L] -Q , 

Jqb>0 



qk'^2nB ^ QbL 



It is of little utility, however, to introduce eigenvalues/eigenvectors of > as the integration 
variables. Rather, it is natural to treat Q^"^ = QbL as a new matrix to integrate over. Properties 
of these (non-Hermitian!) matrices are discussed at length in the Appendix B of [I], and references 

therein. The matrices satisfy 



— LQg ' L, have all eigenvalues real and can be diagonalized 



by a (pseudounitary) similarity transformation: Q^' = TpbT ^, where pB — diag(pi,p2), and 
ub X riB diagonal matrices pi,P2 satisfy: pi > , p2 < 0. Pseudounitary matrices T satisfy: 
f^LT = L and form the group U{nB,nB) ("hyperbolic symmetry"). 

We again introduce the diagonal entries pi and p2 along with the matrices T S u(i)x'^'x'u(i) ^® 
new integration variables. The integration measure dQ^j^^ is given in new variables as |30[ : 

dQ^^' oc dp,dP2 n {p'l'' -p'^'^y {p'^'' -i^'^y n (^^^'^ -p'^-^Ymt) 

where the last factor is the invariant measure on the coset space of T— matrices. 

Again, the only term in the integrand of Eq.js^) which depends on the pseudounitary matrices T 

is obviously the exponential exp A^Tr (^LfiBTpBT^^^ . We immediately see that the corresponding 

integral over the non-compact ("hyperbolic") manifold of T— matrices is exactly that addressed 
by us in Section I (cf. Eq.(||)). 

I{fiB,Pi,P2) = J dt,[f) exp UnTv ( ) ( ) tA (85) 



18 



det 



('l)„('2)l 



(X 



h,l2 = l 



det 



As a final step we change p2 ^ ~P2 and with the integrand depending only on the eigenvalues, 
we arrive to the following expression: 

g-fTrlA.]^ r 

J(^n{P-b, P-f) a r - — 1 A r - — t" / dqpAlqp} {deigp) ^ det 



X e 2 



X det 



2nF 

fei,fc2=l 



-fTr(p?+pi)^g^[^^^^- 



N-2nE 



(il)„(i2) 



2n£ 



det 



2nf 



det [gfel„B - Pi] det [g^l^ 



(86) 



fe=i 



where we denoted i?^^' the integration domain: —oo < qk < oo, k = l,...,2ni? and i?^'' the 
domain < p^'-* < oo for I = 1,2, ...,nB- Taking into account presence of the Vandermonde 
determinant antisymmetric in all q's as well as the symmetry of the rest of the integrand with 
respect to (2np)\ permutations of entries of the matrix qp = diag(gi, 92tif) see that we can 
effectively replace the determinant factor: 

2ni? 



det 



fcl 1^2— 1 



{2npy.e 



and perform similar replacements for the other two determinants in Eq. 

Summing up, we derived the following integral representation for the correlation functions of 
the characteristic polynomials 



ni-^i ZN{-i^i[p)ZN{-iH 2F 



C ' 



dqpAlqp} (detg^) 



2B/ I GUE 



(87) 



A{fiF}A{flB} Jr 

dp,A{p,} dp2A{p2} TT fe'^)+p^'^))e-fTr(p?+p^)^.^Er.\Kip5"-^^'i^^1 



llj2 = l 



2nf 



det [piP2]^ ^""^ '^'^^ i^khiB ~ Pi] det [qkhig + P2] 



k=l 



which is still exact for N >2nB and valid for arbitrary values of parameters such that Im/i^'^ > 
and liJifi^B < 0. 

Before treating the integrals in the limit ^ 00 by the saddle-point method we can restore 
the normalisation constant C by comparing both sides of the equation in the limit N— fixed, 
fi^p"^ 00, Im/Lt^g — > 00, Im/i2^ —00. Obviously, in this limit the presence of the random 
matrix H £ GUE is immaterial and by its very definition the correlation function tends to: 



m=lZNitA%)ZN0^'2L) 



,(0 



(-1) 



AT II 



N 



GUE 



AT 
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(k) 

In the right hand side close inspection shows that the integrals over qp are dominated by vicinity 
of g^'' = iJp ^ . They effectively decouf 
calculated yielding exactly the factor: 



(k) (k) . . (I) . 

of Qp = Hp . They effectively decouple from the integrals over p\ ^ and can be straightforwardly 



A{Af} 



N 



.fe=l 



N 



On the other hand, performing the remaining integrals in the appropriate limit amounts to evalu- 
ating the following expression: 



(-1)" 



dpi 



dp2 A{pi , ~p2}e^ ^'=1 L^^^''^ ' "''^^''^ J det [pip2 



,(0 „(')! 



N-2n 



It can be done by expanding the Vandermonde determinant in the sum over all permutations, 
evaluating the corresponding integrals and resuming the resulting expression back to form another 
Vandermonde determinant: 



-1) 



^b("B+3) 



n 



N-l 

k=N~2ni 



k\ 



N~2nB+l 



Combining all these facts we restore the normalisation constant as: 

(^_l^N{ni+np)-nB{nB/2-l)-j^2nB{N-nB)+nB+n 



(89) 



Coming back to investigating the expression Eq(|^ we can already continue analytically: /i^^* 
i^p "' for k = 1, ...,2np and set all imaginary parts of the spectral parameters /i^'^ and /ijs to zero. 
As usual, we are interested in the so-called "scaling limit" when all the spectral parameters ^J-l^, 
IJ,2B as well as /i^'' are around the same point of the spectrum ^ such that < 2, their mutual 
distance being of the order of N~^. Correspondingly, we set: 



,(0 



(fc) 



N 



and consider all lvb , ujp = 0(1) when N — > oo. 

In this way we reduce the expression under investigation to the form convenient for starting 
the saddle-point analysis: 



Ci 



e 2 



f Tr[Af 



A{lof}A{lub} 



dqi... / dq2nFA{qF} 



'2n F 

X{qk 

.k=l J 



2"F , ,(*:) 



f dp^^i.. f dp^r^ f dp^ii.. f dp("«)A{pi}A{M n fpi'^^+p^' 

Jr^- JRi- Jr+ Jr+ , , ^ 



(90) 



R+ 
'2nf 



X e 



.1=1 



1 9fe - Pi ^ 



k=l 1=1 



qk+P2 
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where 



and 



(91) 
(92) 



Now it is evident that in the Umit N oo the contributions to integrals come from the 
stationary points of the "actions" Cpiq), Cib{p) and Cib{p) given by solutions of the equation 



ijj. 



0: 



(0 
Pi 



*=2 

(0 1 
P2 =2 



k = I, 2np 
I = 1, ...,nB 

Here we took into account the restrictions of the original integration domain: Kep^,^^^. > 



(93) 



Presence of the Vandermonde determinants as well as the factor Jlfcllll Y[i=i 
makes the integrand vanish at the saddle-point sets and thus care should be taken when calculat- 
ing the saddle point contribution to the integral. First of all, the totality of 2^"^ saddle-points 
Qp'' — (qf, q^n^) can be further subdivided into classes giving contributions of different orders 
of magnitude in powers of the small parameter N^^. A little inspection reveals that the leading 
contribution comes from the choice of half of saddle-points to be (j+ , the rest being q~ , with total 
2nF 
np 

brackets inside the Vandermonde determinant Y[i<ki<k2<2nF (^'^i ~ 1k2) is minimal. 

To find the contribution from each of the relevant saddle-point sets explicitly let us subdivide 
the index set 1, 2, 2np into the set {A'+} = (fci < k2 < ... < knp of those indices 1 < fc^ < 2np 
for which qi^^^^ ~ q'^ and the rest denoted as {K^}. Let us also present the integration variables qk 

as : qk£{K±} 



qk - Pi 



(0 



qk 



(0 

P2 



number of such sets 



(compare [|l2[ ^). Indeed, for such a choice the number of vanishing 



"k ' 



from the saddle-point values. Then: 



ith two set of variables — ^a^g^-^^j^ serving to describe deviations 



Mqp} = n ilki~qk2) 

l<ki<k2<2nF 

n (9fci-'7fe2) n iqki-qk^) n i<iki-qk2) (94) 

ki e {K+} fci e {K^} ki e {K+} 

k2 e {K+} k2 e {A'-} k2 e {K^} 

= (-l)'^'+'^'-(4-^2)^A{d+}A{Q!-} + h. o. t. , 

where eK^,K- is odd or even integer serving to take into account the sign factor arising in the process 
of rearranging indices in the last of three products in the above equation. The abbreviation h.o.t. 
stands for higher order terms in a's. Further we expand in the exponentials up to terms quadratic 
with respect to a's and have: 
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k = l '^"'^k 



E.e{K+> E.e{^_} iEfce{x+} ' 



(fc) 



e-^[E::T^-fe)] = (_i)iv„.g-^"-(i+^)-fV^^(«^E.,^._>K)^-«-E.,,.^>K]^^gg^ 

where we made use of q^q~ = — 1 , g"*" + g_ = iji, \ + = ^ ~ and 1 + -j^^ = 

Similarly, we set p^p = q+ + f3^ , p^'^ = — (g" + /3j~) to describe deviations of pfi2) 
their saddle-point values. In this way we obtain: 



2nF riB 



(97) 



fe=l /=! 



[-(4-M^r^"^n n k-at n K-Ai+h.o.t. 

;=i/se{j<:+} fceff+l 
A{p,}A{p2} n (f^i-ii +f2,;.) = (-l)"^ [- (4 - m')] ^ A{^+}A{/3-} + h. o. t.(98) 

ll,l2 = l 



and 



= e 



Er=i EL", -^'ie* Erj^i Ki/'-^+^^'i/T 



e-^Er=i[^-(pi")+^-(f2')] = e-^"B(l+i#)-f^/S^(E,9+[/3^]'-E,9-K]') 



Let us now introduce four diagonal matrices of size n = hb + np- 
e+ = diag(a+,/3+) , 6" = diag (a", /J") 



(99) 



(100) 



O = diag (ci 



0+ = diag (o 



where = diag ( cj^'' : /c € {-ft^±} ) and tjj = diag 



XI) 

hB I ' ""B 



ilig = diag ( cj, 



/ = 1, ...,nB- 



Collecting all the factors we now can represent the leading order contribution to the correlation 
function as A?^ ^ 1 in the form: 



^^(AbjAf) = ^s.p.(As5 Af) 

N 



X J de+A|e+}ox 

X J dQ-A |e-| exp jyg- ^4 - ^u^Tr 



zTr 



X exp 



(-1)-"/X.p.(Ab, Af)(4 - f^r'^A {O-} A 



(101) 



(102) 



2Afx/r 



-Tr 



n- 



-—Tr 
9+ 



0+ 



where 



^s.p.{ij'B,fj-F) oc 



„2 — ATn 



-(-l)^Ci(4-/.^)-e 



(103) 
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„(') 



+ iq 



„(') 



Here we used the integral formula 



deA|e|exp|-^Tr 



e 



iTi 



en 



t— 



'/2 ... r 1 

-A{l^}exp --Tr 



r2 
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and denoted: e = tK+,K^ + Aff^F + nBUp + 

A cl ose i nspection of th e quotient of the Vandermonde determinants occurring when substitut- 
ing Eq.(p^ into Eq.(|Ol]) reveals that: 



A{uf}A{^b} 
where we introduced the notation: 



(fc) , 



[^F -^2b) 



h<l2 [^IB 



OJ. 



2S ) ^ ki£ {K^} 

k2 e {K+} 



— UJp 



(105) 



(106) 



and (— 1)'^^+'^- is exactly the same factor that appeared in our calculation earlier due to rear- 
ranging variables inside the brackets in the product of Vandermonde determinants. 
We also observe that: 



,(fc) 



iq 



(fc) . - 

UJp' -iq 



E 



,(fc) 



fee{A'+} fee{K-} fcefAT-} fce{A:+} 

The last relation allows us to write down the final result of the calculation in the form: 



K^N^ocifj'B, A-f) 



tq ' 



(107) 

, ,(*=) 



+ } 



and the summation goes over all possibilities of subdividing the index set 1,2, ...,2nF into two 
index sets and {K^}. Here CN,nB,nF stands for the overall normalisation constant: 



CN,ns,np = (27r)"s(-l)^"«+"F-"f + — 

and we neglected all the terms of the order of 0{N^^) in the exponential to be consistent with the 
leading order approximation. 

Remembering that the mean spectral density of the GUE eigenvalues in the limit of large N 
is given by the Wigner semicircular law: /o(/i) = ^ ■y/4 — /i^ so that = i itp{ijl) we satisfy 
ourselves that for /i = the derived expression coincides with one announced in and obtained 
by a rather different method. 



23 



6 Conclusions 



In the present paper we have demonstrated that the method suggested in [I] allows one to analyse 
the correlation function containing both positive and negative moments of characteristic polyno- 
mials. This technique combines simultaneous exploitation of the standard Hubbard- Stratonovich 
transformation with an integration theorem (see formula (|8l|)). The latter is a new element as 
compared with [I] which replaces the Ingham-Siegel integration used there. The method leads to a 
compact integral representation (|8^) for the correlation function. To study the asymptotic limit of 
large GUE matrices we needed to expose variables amenable to the saddle point treatment. For this 
purpose we had to derive the integration formula extending the Itzykson-Zuber-Harish-Chandra 
integral to the non-compact Kahler manifold U(ni,n2)/T (see expression (||)). Our derivation is 
based on the Duistermaat-Heckman localization principle. 

Our preliminary considerations show that the outlined procedure works well for other ensembles 
of random matrices, in particular for the chiral GUE and non-Hermitian ensembles. In the latter 
cases it requires a non-compact analogue of the integral formula found by Guhr, Wettig and by 
Jackson, Sener and Verbaarschot [Q. The corresponding calculation will be published elsewhere. 
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Appendix A. Diffusion derivation of the Itzykson-Zuber type 
integral on the pseudo-unitary group U(ni, 



Let us consider a diffusion on matrices that are elements of the Lie algebra u(ni,n2). Any such 
(ni -I- n2) X (ni + na) matrix A satisfies the equation: 

= XA\ (108) 

as it follows from the pseudo-unitarity of the group U(ni, na). The corresponding Laplace operator 
invariant under pseudo-unitary transformations acquires the following form: 

ni+n2 

5] 5V5^'+l/2 i-r^ [dVidReA,f + d'/idlmA,,f] (109) 

i—l l<i<j<ni+n2 

where the symbol aij takes the values when 1 < * < j < rii or ni < i < j < n2 and 1 when 
1 < i < ni < J < n2. Once the Laplacian Da is given, it generates the diffusion on matrices 
satisfying Eq.(105). Such a diffusion is described by the heat equation 

^DAi;{A,t)^dt^P{A,t) (110) 

with the initial condition 

^{A,t^O)^(j){A) (111) 
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The solution of the diffusion problem defined above is represented by an integral over the matrices 
B satisfying the condition Eq.(108): 



(27r<)("i+"2)V2 



(112) 



Let us note that the expression Tr{A — B)'^ is not positive definite for the matrices from the Lie 
algebra u(ni, n2). However one can always choose the initial distribution 0(-B) in a way that ensures 
the existence of the integral Eq.(112). Let us further assume that the initial distribution (t>{B) is 
invariant under pseudo-unitary transformations: 



0(5^5" 



(b{B), ge U(ni,n2) 



(113) 



This condition implies that the function (t){B) depends only on the set eigenvalues of the matrix B 
and is symmetric under separate permutations of the first rii and the rest of n2 eigenvalues. Indeed, 
the corresponding Weyl group is S 

ni ^ and not Sn^+n2 the case of unitary transformations. 
In what follows we adopt the argumentation used in the original work by Itzykson and Zuber 
[p^ to the present case. Let us take a diagonal matrix A in Eq. ( |112| ) and diagonalize the matrix 
B by a pseudo-unitary transformation: 



A = diag(ai,Q;2, • • • ,ani+n2) , B = u diag (/3i, ^2, • ' ' , /3ni+n2) « \ v<EU{ni,n2) 



(114) 



Then the integral expressing the solution of the diffusion problem (Eq.(112)) can be rewritten as 
follows: 



1 



(27ri)("i+"2)V2 



dv / dpA^{f3)exp{~l/2tTT[{a-vPv-^f])<j){P) 



(115) 



where a, (3 stand for the diagonal matrices diag (ai, a2, • • • , ani+n2) and diag {f3i, P2, ■ ■ ■ , /3ni+n2)j re- 
spectively. We introduce the function ^(a, t) antisymmetric with respect to separate permutations 
inside the sets (ai, 02, • • • , a„^) and (apj+i, ani+2, • • • , o:„^+n^): 



C(a, t) ^ A{a)'ilj{a, t), ^{a, t ^ 0) = A{a)(f){a) 



(116) 



Act on this function by the Laplacian operator Dj^. For the function (^(a,t) depending only on 
the eigenvalues of the matrix A and ip{a, t) being a solution of the diffusion problem Eq.(llO), the 
procedure yields the following differential equation: 



ni+r!2 



(117) 



The only diffusion kernel K(a,/3, i) corresponding to the above equation which is antisymmetric 
with respect to separate permutations inside the sets (/3i, /32, • • • , Pm) and (/^m+i, /3ni+2j • • • , /3ni+n2) 
is given by 



(27rt)("i+"2)/2 



PGS„j xS„2 



■A 



P\2 



(118) 



Comparison of the equations Eqs.(|ll5|,|ll6|) and Eq.(118) yields the desired formula Eq.(^) after a 
simple manipulation. 
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Appendix B. Complex parameterization of U(2)/U(l) x U(2) 
and U(l,l)/U(l) X U(l) 

First we note that the manifold U(2)/U(l) x U(l) is equivalent to SU(2)/U(1). The complex Lie 
algebra corresponding to the group SU(2) has three elements e^, e_g, h in its Cartan-Weil basis. In 
the two-dimensional fundamental representation these basis elements are expressed as follows: 



Following the general method of constructing the complex parameterization we decompose the 
representative g{z,z) of the coset space SU(2)/U(1) in the same way as it is done in the formula 
Eq.(^. The factors u{z) and p{z, z) are given by: 

it(2) = exp (zBg) , p(z, z) = exp z)eg) • exp (fc(z, z)/i) (120) 

Any element of a coset of a unitary group must satisfy the unitarity condition. For our case 
the unitarity condition g'^{z,z) = g~^{z,z)yg{z,z) € SU(2)/U(1) is equivalent to the following 
algebraic relation: 

piz,z)pHz,z)^u~\z){uHz)y' (121) 

Thi s rela tion enables one to derive explicitly the functions y{z, z) and fc(z, z) entering the formula 
Eq.(120). We chose to perform the calculations below in the fundamental matrix representation 
since the obtained expressions for y(z, z) and fc(z, z) are the same in any representation. In the 
two-dimensional matrix representation we have: 



«(») = ( J 1 ) , exp(y(z.l)e-,) = ( J 



Inserting the above matrix expressions to the formula ( |121| ) we find: 

i?e(fc(z,z)) =ln(l + zz), y{z,z) ^ z/{l + zz) (123) 

Let us note that the function fc(z, z) is specified up to an arbitrary complex part. This means that 
the corresponding element of the coset space g{z,z) is determined up to a multiplication by a torus 
element from the right, as it must be the case. 

The parameterization of the coset space U(l,l)/U(l) x U(l) is obtained by following similar 
steps. The difference is that the corresponding representative of the coset space should be an 
element of the pseudo-unitary group U(l,l) rather than of U(2). Therefore the representative 
g{z,z) of the coset space U(l,l)/U(l) x U(l) must satisfy the pseudo -uni tary condition Eq.(||). 
The latter leads to the following algebraic relation (which replaces Eq.(121) above): 

p{z, z)\p\z, z) = u-^{z)\{u^)-^ (124) 

Using this formula we obtain the expressions for the real part of the function k{z, z) and for the 
function y(z, z): 

Re{k{z,z)) =\tl{1~ zz), y{z,z) = -z/{l-zz) (125) 

A remarkable feature of the described parameterization is that the expression for Re {k{z, z)) can 
be considered as fundamental Kahler potentials. This is a quite general property common to any 
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homogeneous Kahler manifold with a unitary or pseudo-unitary transformation group as was shown 
by Itoh, Kugo and Kunitomo [Q 

An ahernative way to find the Kahler potentials in terms of the local complex parameters z, z 
is to exploit the relation Eq.(|39|) (or its analogue Eq.(^) for a pseudo-unitary coset). The number 
of fundamental Kahler potentials is equal to the number of projection matrices or to the rank of 
the Lie algebra under consideration. The Lie algebra sm(2) has one basis element h in its Cartan 
subalgebra, so the rank of sm(2) is equal to unity. The projection matrix corresponding to the 
basis element h is determined from the equations (BSh: 




1 



(126) 



Now insert the projection matrix r\ and the matrix u{z) given by Eq.(122) to the formula Eq.(p9[) 
for the unitary coset space U(2)/U(l) x U(l) ( or to its counterpart Eq. ([45| for the pseudo-unitary 
coset U(l,l)/U(l) X U(l)). A simple calculation yields the corresponding Kahler potentials: 

^u(2)/u(i)xu(i)(2^,^) =ln(l + zz), iiru(ia)/u(i)xu(i)(^,^) = ln(l - zz) (127) 

Once the Kahler potentials are known, the (1.1) forms on the manifolds U(2)/U(l) x U(l) and 
U(l,l)/U(l) X U(l) can be immediately obtained from the relation Eq.(||): 



\ dz l\dz \ dz l\dz , , 

f^U(2)/U(l)xU(l)=-^(^^^-)2. f^U(l,l)/U(l)xU(l) = ^(Y3^ (128) 



Finally let us determine the momentum maps on the manifolds U(2)/U(l) x U(l) and U(l, 1)/U(1) x I 
As it can be seen from Eq.(^) the momentum maps are completely determined by the matrix 
p(z, z). In order to find that matrix we insert the projection matrix r\ and the matrix m(z) to 
the formula Eq.(|4l|) for the compact coset U(2)/U(l) x U(l) and to the formula Eq.(|6|) for its 
non-compact counterpart U(l, 1)/U(1) x U(l). We obtain: 

1 / zz z \ — 1 / zz z \ 

PU(2)/U(l)xU(l)(z,^) = ^ PU(l.l)/U(l)xU(l)(^,^) = z -\) ^^^^^ 

Now the momentum maps Tq{z,z),T-q{z,z),Th{z,z) corresponding to the basis elements 6^,6^ 
and h can be easily constructed. For the space U(2)/U(l) x U(l) they are given by 

Tq{z,z) ^ - ^ _ , T_,(z,z) = — Th{z,z) = I , ^! (130) 
1 + zz 1 + zz 1 + zz 

The corresponding momentum maps for the non-compact coset U(l,l)/U(l) x U(l) are 

Appendix C. Transformation of projection matrices pi{z,z) 

In order to prove Eq.(^) we note that the decomposition gu{z) = u{gz)p{z, g) that defines the 
group action on the flag manifold under consideration leads to the following expressions for u{gz) 
and u^{gz): 

u{gz) ^ gu{z)p~'^{z,g), v) {g z) = {p\z , g)) ^ u\z)g'^ (132) 
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Rewrite pi{gz,gz) explicitly using formula Eq.(]4l|): 

Pi{gz,'gz)^u{gz)rji{rjiU^{gz)u{gz)rii + I ~rji) ^ r]iu\gz) (133) 



Insert u{gz) and u''{gz) given by Eg. ( 132 ) to the above formula and use the properties of the 
projection matrices r/j (equation (|3^)). Taking into account relations 

ii~m)ii-m) = ii-m) (134) 

we obtain: 

Piigz.'gz) = 
gu{z)p^^{z,g)r]i (riip^^{z,g)-q^ + I - r],) x 

{r]^u\z)u{z)r|^ + I - iji) ^ (r]^{p^z,g)) ^ + I - rj^ r]i{p\z,g)) ^ u\z)g'^ (135) 



The transformation law for the matrices Pi{z, z) (see Eq.(43)) follows immediately when we simplify 
the above expression using 

P'^{z,g)Vi{VtP^^{z,9)Vt + I - Vi) ^ =Vi, 
{vi{p\z,g)) %^ + /-r/,^ r]^{p^{z,g)) ^ = rj, (136) 

Appendix D. Proof of the Theorem I 

In this Appendix we give a proof of the statement of the Theorem I. In fact, we demonstrate the 
validity of the closely related 
Theorem la 

Consider a function i^(Si, S^) of iV-component real vectors S( 1 <l <m such that 

/ dSi... / d'S™|^^(Si,...,S^)| <oo (137) 

Denoting S"^ the transposition suppose further that the function F depends only on m{m + l)/2 
scalar products S^S;^ 1 < '1,^2 < so that it can be rewritten as a function J-{Qm) of to x m real 
symmetric matrix Qm'- 

Then for > to the integral defined as 

i'-Pm^ I dSi... f dS™F(Si,...,S™) (138) 

is equal to 



'Q,„>0 

where the proportionality constant is given by 



f-ir\ < \ f . / . ^ (Ar-m-l)/2 

W=4m/ rfQmfdetQ™ :F{Qm) (139) 

JQ,„>0 ^ ' 



N,rn j-jm-l p ^_/V^^ 
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and the integration in Eq.( |l39| ) goes over the manifold of real symmetric positive definite m x m 
matrices Q^. 
Proof 

We prove the statement by induction in m for any N > m. 

First, for m = 1 we parameterize S = r (Cjve) , where r — ||S|| > 0, the matrix On € 0{N) is 
a real orthogonal N x N satisfying: OJ^On = Ijv and — 1 first components of the vector e are 
chosen to be zero, the last component being unity: e-^ = (0, 0, 1). The integration measure can 



be written as dS 
0(7V), such that: 



drdii{0), the last factor standing for the Haar's measure on the group 



f T,N/2 



T{N/2) 



N 



Now for iV > 2 we have: 



/ dSF{S')^nNf drr^-'F{r^) = ^ f dqnq[^-'^^' Hln) 
Ji?" Jr+ ^ Jr+ 



which proves the statement and gives the value = as required. 

Suppose now that the statement is true for (m— 1) vectors, each with {N — 1) real components, 
that means: 



^Af-l,m-l ~ -'Af-l.rn-l 



N-l,m-l 



i>0 



dQm-1 (^dctQ 



m — 1 



(N-m-l)/2 



Q 



rn — 1 



(140) 



To consider the case of m vectors, each with N components we represent the m x m matrix 



Si 



Now parameterize = (Oatg) as before, and for fc = l,2,...,TO — 1 introduce new vectors 
Sfc = OjvSfc as integration variables. Obviously, the entries of the matrix Qm-i do not change, 
whereas 



Sk Sm I — 



where 5'jv,fc stands for the last (iV— th) component of the vector S^. Further, let us consider first 
N — 1 components of the vector Sfc as forming the vector (k, for the last A^— th component of the 
vector Sfc using the notation: qk,m., k = 1, ...,m — 1. Obviously: 



dSfc 



dSfc = dCkdqk-. 



Qm— ] 



k.l 



Q 



c 

rn—1 



k,l 



1< fc, / < m - 1 



(141) 
(142) 



where (g^_,)^^^ = (CjO) 

After all those preparations we change the order of integrations (which is legitimate in view of 
the condition Eq.(137) and the Fubini theorem) and represent the integral in the right-hand side 
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of Eq.(|T3|) as: 



In,m^^n I drr^ ^ / rfgi,m-.-dgm-i,m / dCi... / dCm^i Qc_ 

IR+ JR^-^ jRr^'i JRN-I 



where 



Qc 



kA 



Qk,'mQl,m 



J 



Now we can apply the equation Eq. (|l40| ) to replace the integration over the vectors to 
that over the corresponding positive definite matrices. Further introducing quantities qk,m = 
i'm<ikm , fc = 1, TO — f as integration variables and denoting rf-^ = qm.m we immediately see that 
the above integral can be written as: 



N~l rn—l I dqrn,m / dqi „i- 
R.+ JR. 



dQm-iq 



(Af-n-l)/2 



i>0 



detQ,, 



dqra- 

(N-n-l)/2 



(143) 



R 



Qr. 



where we denoted 



qi,r. 

91, r: 



\ 9m, 1 

Using the determinant identity 



(Qm-l)fc^; + ■ 
'?m,2 



dot 



we see that: 



det 



= 5 X det Qrn-l 



X det 



A little more thinking shows that the conditions qm,m > and Qm-i > ensure that Qm is 
a general symmetric positive definite of dimension to. For, all minors of the matrix Qm either 
just coincide with the minors of the matrix Qm-i or with those of the matrix with elements 



-\k,i , qk,mqi,m ^ i^Q^jj being positive definite. 



Combining all this knowledge and the fact that cj^^C^^j^ 



^N\n that the formula 



Eq.(143) can be written exactly as the right-hand side of the equation Eq.(139) thus completing 
the proof. 

Theorem I then follows by trivially repeating, mutatis mutandis, all the steps of the proof given 
above for the case of complex vectors and Hermitian matrices, replacing the orthogonal matrices 
with unitary one as appropriate. 
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